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Abstract 

In this paper, to begin with, we review six different analytical methods which are 
widely used to derive symmetries, integrating factors, multipliers, Darboux polynomi¬ 
als and integrals of second order nonlinear ordinary differential equations. We illustrate 
the algorithm behind each method by considering a nonlinear oscillator equation as an 
example. In the second part of this paper, we examine the interconnections between 
these methods. We establish a road map between extended Prelle-Singer procedure with 
all other methods cited above and thereby demonstrate the interplay between Lie point 
symmetries, A-symmetries, adjoint symmetries, null forms, integrating factors, Darboux 
polynomials, and Jacobi multipliers of second order integrable ODEs. The interconnec¬ 
tions are illustrated with the same example finally. 

1 Introduction 

During the past few years there has been a great deal of interest in finding the interrelation 
between symmetries, integrating factors, integrals of motion and integrability of the associated 
dynamical systems, especially for finite degrees of freedom systems [TJ. Confining our attention 
on dynamical systems described by second order linear/nonlinear ordinary differential equa¬ 
tions, the integrability is guaranteed either by providing a time independent integral (which 
can be interpreted as the Hamiltonian) or two time dependent integrals (2j. The required 
integrals can be constructed with the help of symmetries (it may be Lie point symmetries, 
adjoint symmetries, A-synnnetries, contact symmetries or nonlocal symmetries), integrating 
factors and multipliers. To investigate these quantities, that is symmetries, integrating factors 
and integrals, several powerful methods have been developed in the literature over a long pe¬ 
riod 1312 El- The methods that are widely used to explore these quantities are the following 
(chronological order): (i) Jacobi last multiplier, (ii) Darboux polynomials, (iii) Lie point sym¬ 
metries, (iv) adjoint symmetries, (v) A-symmetries and (vi) Prelie-Singer method. We note 
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here that we have listed out only the popular methods in the literature and left the others 
which are applicable only for specific equations. 

Since we are going to concentrate only on these methods, in the following, we recall briefly 
the development of each method. The Jacobi last multiplier method was introduced by Jacobi 
in the year 1844 mm- The method lay dormant for several decades. Only recently Nucci 
and her collaborators have demonstrated the applicability of this method in exploring non¬ 
standard Lagrangians associated with certain second order nonlinear ODEs EElina. One of 
the important applications of multipliers is that one can determine the integrals associated 
with the given equation by evaluating their ratios. 

Darboux polynomials method was introduced by Darboux in the year 1878 DU • There exist 
a lot of different names in the literature for Darboux polynomials, for example we can find: 
special integrals, eigenpolynomials, algebraic invariant curves, particular algebraic solutions 
or special polynomials. The method also provides a strategy to find first integrals. If we 
have enough Darboux polynomials then there exists a rational first integral. Darboux showed 
that if we have d ^- d ^ 4 - 2 Darboux polynomials, where d is the order of the given equation, 
then there exists a rational first integral which can be expressed in terms of these polynomials 

mmmm- 

Lie symmetry analysis is one of the powerful methods to investigate the group invariance 
properties of the given dynamical system, and establishing the integrability of the given dif¬ 
ferential equation is just a by-product of this method. The method was developed by Sophus 
Lie at the end of nineteenth century. In this method one first identifies the symmetry vector 
fields associated with the differential equation. The symmetry vector fields can then be used 
to derive integrating factors, conserved quantities, similarity variables and so on [3]. 

Exploring integrating factors and obtaining integrals for the given equation is one of the 
classical ideas and there exists several methods in the literature to obtain the integrating 
factors. One such method is the adjoint symmetry method. This method was developed by 
Bluman and Anco in the year 1998 [5]. 

All the nonlinear ODEs do not necessarily admit Lie point symmetries. Under such a 
circumstance one may look for generalized symmetries (other than Lie point symmetries) 
associated with the given equation. One such generalized symmetry is the A-symmetry. The 
A-symmetries can be derived by a well defined algorithm which includes Lie point symmetries 
as a specific subcase and have an associated order reduction procedure which is similar to the 
classical Lie method of reduction [5]. Although A-symmetries are not Lie point symmetries, 
the unique prolongation of vector fields to the space of variables (t, x ,..., x n ) for which Lie 
reduction method applies is always a A-prolongation for some function A (t,x,x). For more 
details, one may see the works of Muriel and Romero [4j. 

Recently, a powerful method, namely the extended Prelie-Singer method was developed 
by our group jl 6 | . The method is an extended and complete picture of the technique given 
by Duarte et al. for second order ODEs nu. The original method, which is applicable to 
planar ODEs, was proposed by Prelie-Singer in the year 1983 [18| . The attractiveness of 
the Prelie-Singer (PS) method is that the method guarantees that a solution will be found 
if the given system of first-order ODEs has a solution in terms of elementary functions. The 
same methodology has also been demonstrated in second and higher order nonlinear ODEs 
including coupled ones [ 16] . 

Even though the aforementioned methods have been introduced at different times in the lit¬ 
erature recently attempts have been made to investigate the interconnection between symme¬ 
tries, multipliers, integrating factors and integrals as well as the methods that produce them. 


2 



Upon analyzing the factors the interconnection have been found between (i) A-symmetries 
and Lie point symmetries [3J, (ii) Jacobi last multiplier with Darboux polynomials P2L (iii) 
A-symmetries and null forms |4j and (iv) Lie point symmetries and Jacobi last multiplier [8j. 
A natural question which arises here is whether there exists a more encompassing intercon¬ 
nection which relates all these methods/quantities. In this paper we establish a road map 
between extended Prelle-Singer procedure with all other methods cited above and thereby 
demonstrate the interplay between Lie point symmetries, A-symmetries, adjoint symmetries, 
null forms, integrating factors, Darboux polynomials, and Jacobi multipliers of second order 
integrable ODEs. For more details, see our recent work |20| . 

We organize the presentation as follows: In Section 2, we give an overview about these 
methods with an example. In Section 3, we give the interlinks between these methods. We 
also demonstrate the interlinks with the same example. Finally, we present our conclusions in 
Sec.4. 

2 Analytical methods 

In this section we briefly recall six analytical methods which are widely used in the contem¬ 
porary literature to derive the following quantities, namely (i) integrals, (ii) symmetries, (iii) 
integrating factors and (iv) multipliers. To make our presentation concise and comparable we 
consider the same example in all the methods and illustrate the ideas behind each method. 

2.1 Jacobi last multiplier (JLM) |6j, [TJ 

Let us consider a second order ODE 


x = x, x), 


(1) 


where <f> is an analytic function of the variables t, x and x. Let us rewrite the second order 
ODE (H|) into a system of two first-order ODEs: 

Tt =!{x ' y) ' Tt =3 ^ v) - <2) 

Dividing the second expression by the first expression of Q and multiplying the resultant 
equation by an integrating factor M(x,y ) we arrive at 


M(gdx — fdy) = 0 =>- dl = 0. 


(3) 


Let us assume that Eq.( 
this case, we have 


has a first integral I(x,y ) which is constant on the solutions. In 
dl dl 

dl = —dx + —dy = 0. (4) 

ox oy 


Comparing Eqs. (J4|) and ([3]), we find ^ = Mg and = —Mf. By imposing the compatibility 

condition, q x q v = g y g x , we get M y g + Mg y = —M x f — Mf x . The later equation can be 
rewritten as 


D[M] 

M 


+ fx + 9y 


0 => D[log M] + fx + g y = 0, 


^ d d d 
D = ¥t +x fo +y d^- 


(5) 
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Choosing the functions / and g suitably, the determining equation for the Jacobi last 
multiplier for a second order ODE can be written as 

D[logM\+cj) i = 0. ( 6 ) 

Solving Eq.© we can obtain the Jacobi last multipliers M. Once sufficient number of Jacobi 
last multipliers are determined the associated first integral can be found from their ratios. 

To demonstrate that the ratio of two multipliers defines an integral, let us assume that 
M\ and M 2 are two Jacobi last multipliers. In this case, we have 

ZJ[log Mi] + = 0, D[log M 2 ] + = 0. (7) 

Let I be the ratio of two Jacobi last multipliers, that is I = -jgC Evaluating the total derivative 
we find 

dI_dfM x \_ _ M ± M 2 - M X M 2 _ Mi M 2 
dt~dt\W 2 ) ~ “ Mf ~Wx~W 2 

= —— log Ml - 4; log M 2 = (px - (px = 0. 

at at 

Thus if we have sufficient number of multipliers we can obtain the necessary integrals. 


2.1.1 Example 


Let us illustrate the above method with the following second order nonlinear ODE, namely 

q r 2 

x - : -b 2x 3 = 0 . ( 8 ) 

2 x 

The Jacobi multipliers of Eq.© can be determined by integrating the following equation, that 
is M t + xM x + (|| — 2 x 3 )M± + ^M = 0 (vide Eq.©). To obtain particular solutions of this 
equation we assume the following ansatz for M, namely M = a ( t x )+b(t,x)±+c(t x)±‘ 1 > w l iere 
and c are arbitrary functions of their arguments. Substituting M and its derivatives in the 
determining equation for M and rearranging the resultant expression we obtain a polynomial 
equation in x. Equating the various powers of x to zero, we obtain a set of differential 
equations for the unknowns a, b and c. Solving them we can find the explicit forms of a, b and 
c. Substituting them back in the expression for M we can get 


Mi 


M 3 


x 2 + Ax 4 ' M 2 2xx + t(x 2 + Ax 4 ) ’ 

-Tj-, M 4 = —T. 

2txx + 2x 2 + i 2 (%- + 2x 4 ) x 


(9) 

( 10 ) 


The ratio between the JLMs Mi and M 4 gives an integral (Ji) and the ratio between the 
multipliers M 2 and M 4 yields another integral ( I 2 ) of the form 


h 


x 2 

—7 + Ax, I 2 
x J 


Axx + 2t{Ax 4 + x 2 ) 


( 11 ) 


One can straightforwardly check that = 0, i = 1,2, and the integrals I\ and I 2 are 
functionally independent. From these two integrals we can derive the general solution of © 
in the form 


I 6/1 

64 + if - Ahl 2 t + 4/ft 2 ' 


( 12 ) 
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We note here that all other ratios of M[s , i = 1,2, 3,4, provide only functionally dependent 
integrals, that is the resultant integrals turn out to be either a function of I\ or I 2 or their 
combinations, and so no further information is obtained. 

2.2 Darboux polynomials [llj 

Let us consider a second order ODE (]T]) which admits a first integral of the form I = 
where F and G are functions of their arguments. Taking a total derivative and rewriting the 
resultant expression, we get 

^ ^ = 0 => F = g(t , x, x)F =>• D[F] = g(t, x, x)F, (13) 

where D is the total differential operator and g(t, x,x) = ^ is the cofactor. The above 
equation is the determining equation for the Darboux polynomials nu. Solving Eg. (1131) . we 
can obtain Darboux polynomials ( F ) and the cofactors ( g ). Once sufficient number of Darboux 
polynomials are known the associated first integral can be determined as follows. The ratio of 
two Darboux polynomials with the same cofactor defines the first integral HU. For example, 
let /1 and f 2 be two Darboux polynomials with the same cofactor, say g(t,x,x). In this case, 
we have 

#]_ _ , dh _ t nA \ 

dt 9 U dt 9 2 ' ^ 

Let I be the ratio of two Darboux polynomials, that is I = jq. From this expression, we find 

*L = ±(Ii\ = n A hh 

dt dt\f 2 ) f 2 f 2 

0 = 9 p-- 9 ^p- = 0. (15) 

J 2 / 2 

Hence in this method also the necessary integrals can be deduced with sufficient number of 
polynomials. The combinations of Darboux polynomials also define a first integral. 

2.2.1 Example 

To illustrate the Darboux polynomials approach we again consider the same example given 
above, that is Eq. (|8]| . The Darboux polynomials can be determined by solving the equation 
Ft + xF x + (§§ — 2 x 3 )F x = gF (vide Ea. (1131) 1. To solve this equation we assume an ansatz 
for F as F = a(t,x) + b(t,x)x + c(t,x)x 2 , where a,b and c are arbitrary functions of their 
arguments. Substituting this form in (1131) and rearranging the resultant equation one obtains 
a polynomial equation in x. Equating the various powers of x to zero we get a set of equations 
for the unknowns a, b and c. Solving them consistently we can find the exact forms of a, b and 
c. Substituting them back in F we obtain Darboux polynomials of (j8j) in the form 


F\ 

= ^(x 2 + 4x 4 ), F 2 = 2xx + t(x 2 + 4x 4 ), 

(16) 

F 3 

. 2 

= 2txx + 2x 2 + t 2 (—+ 2x 4 ), F A = x 3 . 

(17) 


Interestingly all the above Darboux polynomials share the same cofactor (g), that is In 
other words, one can verify that D[Fj] = ^F t , i = 1,2, 3,4. Since the ratio of two Darboux 
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polynomials which have the same cofactor defines a first integral, we can extract the required 
integrals from the above Darboux polynomials. For example, the ratio of F\ and F 4 provides 
an integral which exactly coincides with I\ given in dill) and the ratio between F 2 and F 4 
gives another integral which in turn matches with the second integral given in Eg. dill) . The 
other ratios provide only functionally dependent integrals. From these two integrals one can 
get the general solution (fl2l) . 

2.3 Lie point symmetries [3j 

Let us consider a second order ODE of the form ©• The invariance of Eq.fll]) under an one 
parameter group of Lie point symmetries, corresponding to infinitesimal transformations, 

T = t + e £(t, x), X = x + er](t, x), e<C 1, (18) 

where £(f, x) and rj(t, x) are the infinitesimal point symmetries and e is a small parameter, is 
given by 

i^ + + ^ + ^ x ~ ^ ~~ (% + ( 2 Vtx ~ itt)x + (Vxx ~ 2 itx)x 2 

- CxxX 3 + ( rj x - ‘lit - 3x£ x )x) = 0. (19) 

Substituting the known expression cj> in d!9D and solving the resultant equation we can get the 
Lie point symmetries associated with the given ODE. The associated vector field is given by 

One may also introduce a characteristics Q = ij — and rewrite the invariance condition 
m in terms of a single variable Q as 

d 2 Q , dQ _ 

dt 2 4>± dt <t>xQ- 0 . (20) 

Solving this equation one can get Q. From Q one can recover the infinitesimals £ and ij. 


2.3.1 Example 

To demonstrate the Lie point symmetries approach we again consider the same example given 
in Eq.(jSJ). The invariance of Eq.jH]) under the one parameter Lie group of infinitesimal trans¬ 
formations m is given by 


3i; 


r?(^2 + 6x )~( r lt + A (Vx - it) - X ix)) 


2x 


T Vtt T (2 i~)tx itt)x T (jjxx 2 itx)x 


3 x 


• 2 


-ixxX + (rj x - 2 it - 3 xi x )(- - 2x 6 ) = 0. (21) 


2x 


Equating the coefficients of various powers of x m ,m = 0,1,2,3, to zero and solving the 
resultant partial differential equations we can find the explicit forms of i and ij which turn out 
to be i = b\ + + b^t 2 and ij = —( 62 ^ + 2bstx), where b(s, i = 1,2, 3, are arbitrary constants. 

The associated vector fields can be written as (V = ij ^ + r im) 


Vi = 


d t r d d o 9 n d 

at' V2 = t m~ x ai’ V3 = t ai- 2tx T x 
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2.4 Adjoint symmetries [5] 

In general, for systems of one or more ODEs, an integrating factor is a set of functions, 
multiplying each of the ODEs, which yields a first integral. If the system is self-adjoint, then 
its integrating factors are necessarily solutions of its linearized system. Such solutions are the 
symmetries of the given system of ODEs. If a given system of ODEs is not self-adjoint, then 
its integrating factors are necessarily solutions of the adjoint system of its linearized system. 
Such solutions are known as adjoint symmetries of the given system of ODEs [5]- The adjoint 
ODE of linearized symmetry condition (12U1) can be written as 

^ + |(^A)-0 X A = O. (23) 

Evaluating Ea. (l23|) and collecting the powers of x, we get 

A tx T A^x -|- 2Aa; -f- A.(f>xx T ‘^4 > x Ax T (f>A- xx = 0) (24) 

A tt T 2A tx& T k xx x T A.cj)fx T A (f>x X x T A 4^x^-xX 

A(^> 3 , A x4^ + T (fr^-xxX T A x (pt T A x (p x x — 0. (25) 


Solutions of the Ea. (|25|) are called adjoint symmetries. If these solutions also satisfy Ea. (|24|) 
then they become integrating factors. Once we know the integrating factors, we can find the 
first integrals. Multiplying the given equation by these integrating factors, and rewriting the 
resultant equation as a perfect differentiable function, we can identify the first integrals, that 

is 


A i(t, x , x)(.x 


x , x)) = —Ii, i — 1,2, ...m. 


(26) 


2.4.1 Example 


If the symmetry determining equation (1201) is the same as its adjoint determining equation 
(1231) . then the adjoint-symmetries are symmetries. In this case, the equation is called self- 
adjoint equation. If the given equation is not self adjoint, then the solutions of Eo. (l23() are 
called adjoint symmetries. Since the example ((21) is not self-adjoint we have to find the adjoint 
symmetries. So we have to solve the following differential equation + ^(^A) + (j% + 
6j; 2 )A = 0. To achieve this task we assume an ansatz for A as A = h ^ t x)x i ' 

where a[s and b'-s,i = 1, 2, j = 1,2,3, are all to be determined. Substituting this ansatz in 
(1231) with 4> = — 2x 3 and solving the resultant equation we can get the following particular 

solutions, that is 


Ai = 


—2x 


A2 


— (x + tx) 

4 tx 4 + 2xx + tx 2 ' 


(27) 


We find that these two expressions also satisfy Ea. (|24|) and so they become the integrating 
factors of (J8j) as well. Multiplying the given Eq.([8|) by these integrating factors separately and 
rewriting them as a perfect derivative we can find the same integrals I\ and I 2 given in dill) . 


2.5 A-symmetries approach [4] 

All the nonlinear ODEs do not necessarily admit Lie point symmetries. Under such a cir¬ 
cumstance one may look for generalized symmetries (other than Lie point symmetries) asso¬ 
ciated with the given equation. One such generalized symmetry is the A-symmetry. Let V = 
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£(i, x)-^+r](t, x) be a A-symmetry of the given ODE for some function A = A (t, x, x). The in¬ 
variance of the given ODE under A-symmetry vector Held is given by yDA 2 )] (£ —^(t, x, x)) = 0, 
where I/t A A 2 )] is given by + r/-g | + r/^A 1 )]^ 4 - ?y[ A A 2 )]^. Here r/^A 1 )] and r/ A A 2 )] are fi rs t 

and second A- prolongations. They are given by 

h^A 1 )! = (D t + X)r] [x ’^(t,x) - (D t + A)(£(t,x))x, 

ryAA 2 )] = + A) 7 y [A ’ ( 1 ) 1 (f,x,x) - (D t + \)(£(t,x))x. (28) 

Expanding the invariance condition, we find 

£<l>t + ri&x + - ?/ A ’ (2)1 = 0. (29) 


If we put A = 0, we get the Lie prolongation formula. Solving the invariance condition (I29|) 
we can obtain the explicit forms of £, 7 ] and A. Once the functions £, r] and A are known the 
integrals can be evaluated using the following procedure [3]. 

1. Find a first integral w(t,x,x ) of nl A A 1 )] ) that is a particular solution of the equation 

w x + \w± = 0, (30) 

where subscripts denote partial derivative with respect to that variable and n^A 1 )] i s the 
first order A-prolongation of the vector field v. 

2 . Evaluate A(w) and express A{w) in terms of (t, w) as A(w) = F(t,w). 

3. Find a first integral G of dt + F(t, w)d w . 

4. Evaluate I(t,x,x ) = G(t,w(t,x,x)). 


If the given ODE admits Lie point symmetries, then the A-symmetries can be derived without 
solving the A-prolongation condition. In this case the A-symmetries can be uniquely deduced 
from 


_ D[Q\ 
Q 


(31) 


where D is the total differential operator and 

by v = 


Q = 7 / — £x. The associated vector field is given 


2.5.1 Example 

To illustrate the method we again consider the same example, namely Eq. (| 8 ]). In the first step 
we determine the A-symmetries of this equation. Since we already know the system admits 
Lie point symmetries, one can capture the A-symmetries straightforwardly by recalling the 
relation (13111 . The resultant expressions read 

_ 3x 2x 3 _ 4xx + 3 tx 2 — 4 tx 4 _ 1 3x 2 tx 3 . . 

1 2 x x ’ 2 2 x 2 + txx ’ 3 t 2x 2x + tx' 

With £ = 0 and 77 = 1, one can easily verify that the invariance condition (1291) is satisfied for 
all the above three forms Ai, A 2 and A 3 . The associated vector field is given by V = In the 
second step we have to derive the integral associated with each A-symmetry which can be done 





through the algorithm given above. For illustration purpose we consider only one case. For 
the remaining two cases one can adopt the same procedure and derive the associated integrals. 


Substituting Ai = || — 2|i in (13U1) one gets 


w x + 


3x 

2 x 


2 x 3 

x 


Wx = 0 . 


(33) 


To obtain a particular solution of (133[) we solve the associated characteristic equation, namely 


dx 


dx 


\ 3x _ 2.r’ ' 

2x x 


(34) 


Integrating (1341) we find a first integral w(t,x,x ) which is of the form 


w(t, x,x) = —t + Ax. 


(35) 


In the second step one has to evaluate A(w), that is 
. . , dw dw ,3x i.dw 

A( ” ) = W + I & +( 2 


A(w) = F = 0. 


(36) 


Since F = 0, the function w becomes the integral I\. In other words we have I± = ^ + Ax. 
The integral I\ exactly matches with the one given in Eg. dill) . In the same way, one can get 
the other two integrals from the A 2 and A 3 , respectively of the form 


I2 = 


Atx 4 + 2xx + tx 2 




h = 


Ax 2 + At 2 x 4 + Atxx + t 2 x 2 


(37) 


Here also one can check that the total derivatives vanish, that is = 0, i = 1,2, 3. In the 

above I\ and I 2 are the two independent integrals and I 3 = ~ 2 ^—- The general solution (I12|) 
can be found from the integrals I\ and I 2 given above. 


2.6 Prelle-Singer method [16, [T7, IT8] 

Let us assume that the ODE (|TJ) admits a first integral I(t,x,x ) = C, with C constant on the 
solutions, so that the total differential gives 

dl = Itdt + I x dx + I±dx = 0, (38) 

where the subscript denotes partial differentiation with respect to that variable. Rewriting 
equation ([T]) in the form cfrdt — dx = 0 and adding a null term S(t,x,x)x dt — S(t,x,x)dx to 
the latter, we obtain that on the solutions the 1 -form 

dt — Sdx — dx = 0. (39) 

Hence, on the solutions, the 1-forms (l38l) and (1391) must be proportional. Multiplying (1391) 
by the factor R(t,x,x ) which acts as the integrating factor for (1391) . we have on the solutions 
that 



dl = R{(j> + Sx)dt — RSdx — Rdx = 0. 


(40) 
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Comparing equations (1381) with (l40l) we end up with the three relations which relates the 
integral(J), integrating factor (R) and the null term(S'), 


It = R(<t> + xS), I x = -RS, h = -R. (41) 

The compatibility conditions between these relations that is It x = Ixt , It± = Ixx = I xx 
lead us to the following determining equations for the unknowns R and S, namely 

-D[S] = —i i> x + S<j>x + S 2 , (42) 

D[R}= -R(S + </>i), (43) 

R x = R±S + RS±, (44) 


where D = + x-^ + Equations (|42|> - (144[) can be solved in principle in the following 

way. From ([42 [I we can find S. Once S is known then (|43l) becomes the determining equation 
for the function R. Solving the latter, one can get an explicit form for R. Now the functions 
R and S have to satisfy an extra constraint, that is equation (1441) . Once a compatible solution 
satisfying all the three equations have been found then the functions R and S fix the integral 
of motion I(t,x,x ) by the relation 


I(t,x,x) = j R{cf) + xS)dt — 
R+4r\ f R{(t> + ±S)dt - [ [RS + 


dx 


RS + 

— f 

dx I 


d 


dx 

R(4> + xS)dt ) dx 


R{4> + xS)dt ) dx 
dx. 


(45) 


Equation (1451) can be derived straightforwardly by integrating the three relations which relates 
the integral(I), integrating factor (R) and the null term(S'). Note that for every independent 
set (S, R), equation ([45 p defines an integral. Once we find sufficient number of independent 
integrals (for example two independent integrals for a second order nonlinear ODE) we can 
write the general solution. In the following, we illustrate the method with the example (l8|) . 


2.6.1 Example 

To obtain the integrals and the general solution of (]5[), in the PS method, we have to find null 
forms and their associated integrating factors. To obtain the null form we have to solve the 
following equation, namely 


St +xS x + (--- 2x 3 )S x - + S— + S 2 . (46) 

To solve this equation we assume an ansatz for S in the form S = ^p^)+& 2 (t’,xj±+b 3 (t’x)± 2 ’ 
where a(s, and b's , i = 1, 2, 3 are functions of their arguments. Substituting this ansatz in the 
S determining equation and rearranging the resultant equations we get a polynomial equation 
in x. Equating the coefficients of various powers of x, to zero we get a set of partial differential 
equations for the unknown functions a[s and 6(s. Solving them we obtain their explicit forms. 
In our case we find the following two particular solutions, namely 


Si 


2 x 3 

x 


3x 

2 x' 


S 2 


Itx 4 — Ixx — 3 tx 2 
2 x 2 + txx 


(47) 
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Once we know the null forms, substituting them separately in (1431) and solving the resultant 
equation we can get the integrating factors. In both the cases we assume an ansatz for R in 
the form 


ai(t, x ) + a 2 (t, x)x 


7-> _ _ I \ L 1 ) ~ _£ V 2 / _ /^n\ 

b\(t,x) + b 2 (t,x)x + b 3 (t,x)x 2 ’ 

where a[s and bjS, i = 1,2, j = 1,2,3, are arbitrary functions of their arguments. With this 
ansatz we solve the Eq. (I43|) and obtain the integrating factors in the following form 


Rt = R 2 = —Y±F). (49) 

x 3 4fx 4 + 2xx + tx 2 

The expressions (Si, Ri), i = 1,2, satisfy the third equation (I44p also. With the help of these 
null forms and integrating factors, we can find the integrals by substituting them in (1451) and 
evaluating the integrals. Doing so we find the same integrals as given in Eg. (HID . From these 
two integrals one can obtain the general solution (I12D . 


3 Interconnections 

In the previous section we illustrated certain methods which are used to derive integrating 
factors, symmetries and integrals of a given dynamical system. In this section we demonstrate 
from null forms and integrating factors of PS procedure one can recover Lie point symmetries, 
A-symmetries, adjoint symmetries, Darboux polynomials and Jacobi multipliers. To demon¬ 
strate that the null form S and the integrating factor R are intimately related with other 
measures of integrability we identify the following relations. By introducing a transformation 

S = -D[X]/X, (50) 

Ea. (l42p becomes a linear equation in the new variable X, that is 

D 2 [X] = (/> ± D[X] + <I> x X, (51) 

where D is the total differential operator. With another change of variable 

R = X/F , (52) 

where F(t,x,x ) is a function to be determined, we can rewrite Ea. (l43D in a compact form in 
the new variable F as 

D[F\ = 4>±F. (53) 

In the following, we demonstrate that the functions X and F are intimately related to Lie 
symmetries, A-symmetries, Darboux polynomials, Jacobi last multiplier and adjoint symme¬ 
tries. By doing this we also capture the known connections. For more details, see our recent 
work [2Q]. 

3.1 Connection between A- symmetries and null forms 

If we replace S = —Y in equation (142 D we get 

D[Y] = <t> x + Y<j>i-Y 2 , (54) 

which is nothing but the determining equation for the A-symmetries for a second order ODE 
[5] which in turn establishes the connection between A-symmetries and null forms. It has 
also been shown |3j that once Lie point symmetries are known, the A-symmetries can be 
constructed through the relation Y = D[Q\/Q which is also confirmed here. 
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3.2 Connection between JLM and Lie point symmetries 


The Jacobi last multiplier can also be calculated from the Lie point symmetries by evaluating 
the determinant Eld 


A = 


1 

c 

£2 


X 

X 


(1) 

m 

Vi 


(i) 

m 

V 2 


(55) 


where (£ 1 , 771 ) and (£ 2 , V 2 ) are the infinitesimal Lie point symmetries and r/j 1 ' 1 and are their 
corresponding first prolongations. Inverting the determinant we get Jacobi last multiplier 
M = ( 2 ). The determinant establishes the connection between multiplier and Lie point 
symmetries Ed. 


3.3 Connection between Lie point symmetries and null forms 

We have seen that the invariance condition v' 2 ' [x — (f)(t, x, £)] = 0 determines the infinitesimal 
symmetries £ and ij explicitly |3j. The invariance condition in terms of the evolutionary vector 
field Q is given by 

D 2 [Q] = feD[Q] + <l> x Q. (56) 

Comparing Eas. (I51f) and (l56l) we find that 

X = Q. (57) 

In other words the 5-determining equation (15ID now becomes exactly the determining equation 
for the Lie point symmetries £ and 77 with Q = 77 — i;£. Since 5 = — D ^' , the null form 5 can 
also be determined once £ and 77 are known. This establishes the connection between the null 
forms 5 with Lie point symmetries (£ and 77). 


3.4 Darboux polynomials and integrating factors 

Let us consider the function G = where f[ s are Darboux polynomials and n(s are 

rational numbers. If we can identify a sufficient number of Darboux polynomials (irreducible 
polynomials) f[ s, satisfying the relations D[fi]/fi = cti, where a^s are the co-factors, then 

D[G\/G = Yni D[ ! l] = rijai. (58) 

i U 

Now we compare Ea. (l58p with Ea. (l53l) . If we choose G = F then equation (1531) becomes 

D\fi\ a 

2_^ n i—f— = ( Px- (59) 

i h 

In other words the Darboux polynomials constitute the solution of Ea. (l53l) . Since X is al¬ 
ready known, once Darboux polynomials are known the integrating factors can be fixed (vide 
Eq.(E2D). 
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3.5 JLM and Darboux polynomials/integrating factors 

Comparing the equation (l53l) with (J 6 J) we find that 

F = M~\ (60) 

Thus from the knowledge of the multiplier Ad = (-^) we can also fix the explicit form of F 
which appears in the denominator of integrating factor R (vide Eg. (1521) 1 as F = A or vice 
versa. 

3.6 Adjoint symmetries and integrating factors 

Let us rewrite the coupled equations (1421) and (1431) into an equation for the single function R. 
Then the resultant equation turns out to be of the form 

D 2 [R} + D[<f> A R] - 4> X R = 0. (61) 

Comparing the above two equations (1231) and (161 p one can conclude that the integrating 
factor R is nothing but the adjoint symmetry A, that is 

R = A. (62) 

Thus the integrating factor turns out to be the adjoint symmetry of the given second order 

ODE. 

3.7 Example 

We illustrate the interconnections by considering the same example (l 8 |). In the above, we have 
shown that the interconnections between the several analytical methods to solve second order 
differential equations. Since all the quantities are interrelated we start with the Lie point 
symmetries. From the Lie point symmetries, V\ and V2, we can determine a multiplier by 
evaluating the determinant 


1 

X 

X 


1 

X 

X 

= + 4x 4 ). 



Vi 

V? 

= 

1 

0 

0 

(63) 

6 

72 

V2 ] 


t —x 

— 2x 

z 



The reciprocal of this determinant (^A-) gives a last multiplier which coincides exactly with 
Adi (vide Eq.Q). From the Lie point symmetries Vj and V 3 , we can get A 13 and from V 2 
and V 3 we can get A 23 . The expressions and defines the multipliers M 2 and M 3 
respectively which in turn exactly matches with expressions given in Q and COD. From the 
null forms (I47p . we can fix the A-symmetries as S = —Y. The adjoint symmetries are nothing 
but the integrating factors has already been observed through the expressions (I62|) and (1491) . 
The reciprocals of multipliers provide the Darboux polynomials which can be straightforwardly 
verified in this example. 
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4 Conclusion 


In this paper, we have considered six different analytical methods which are widely used in 
the contemporary literature. In particular, we have considered (i) Jacobi last multiplier, (ii) 
Darboux polynomials, (iii) Lie point symmetries, (iv) adjoint symmetries, (v) A-synnnetries 
and (vi) Prelle-Singer method. To begin with we have recalled each one of the methods 
separately with an example. The same example is considered for all the methods in order to 
illustrate the merits and demerits of each method compared to other methods. These methods 
essentially seek either one or more of the following quantities, namely (i) symmetries, (ii) 
integrating factors, (iii) multipliers, (iv) Darboux polynomials and (v) integrals. We have 
investigated the interlinks between these methods in detail. After performing a thorough 
analysis we have shown that starting from the null form and the integrating factor in the Prelle- 
Singer procedure one can obtain A-synunetries, Lie point symmetries, Darboux polynomials 
and adjoint symmetries. We have explicitly shown the interconnections between these factors. 
The same example has been considered to illustrate the interlinks also. 
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